t3 



(N 



CO 

in 



X 



A Peculiar Jet and Arc of Molecular Gas Toward the 
Rich and Young Stellar Cluster Westerlund 2 and a TeV 

Gamma Ray Source 



Yasuo FuKUi/ Naoko Furukawa/ Thomas M. Dame,^ Joanne R. Dawson/-^ Hiroaki 

On " 

O ■ Yamamoto/ Gavin P. Rowell, Felix Aharonian,^''' Werner Hofmann,'' Emma de Oiia 

o ■ 

! Onishi/ Akira MizuNO,^ Shigehiro Nagataki 

^ ■ ^Department of Astrophysics, Nagoya University, Furocho, Chikusa-ku, Nagoya, Aichi, 464-8602 
^Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA. 
^Australia Telescope National Facility, CSIRO, P.O. Box 76, Epping NSW 1710, Australia 
School of Chemistry and Physics, University of Adekaide, Adelaide 5005, Australia. 
^Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland. 
Q^. M ax- Planck- Institut fur Kernphysik, PC Box 103980, 69029 Heidelberg, Germany. 

Q . '^Department of Physics, Faculty of Science, Hokkaido University, N10W8, Kita-ku, Sapporo 



060-0810 

^National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588, Japan. 
^ Solar- Terrestrial Environment Laboratory, Nagoya University, Furocho, Chikusa-ku, Nagoya 



> : 464-8601 
^ \ Yukawa Institute for Theoretical Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502 



fukui@a. phys.nagoya-u.ac.jp 



O ' (Received 2000 December 31; accepted 2001 January 1) 

O ■ Abstract 



We have discovered remarkable jet- and arc-like molecular features toward the rich 
and young stellar cluster Westerlund2. The jet has a length of ~ 100 pc and a width 
of ^ 10 pc, while the arc shows a crescent shape with a radius of ^30 pc. These 
molecular features each have masses of ~ 10^ Mq and show spatial correlations with 
the surrounding lower density HI gas. The jet also shows an intriguing positional 
alignment with the core of the TeV gamma ray source HESS J1023-575 and with the 
MeV/GeV gamma-ray source recently reported by the Fermi collaboration. We argue 
that the jet and arc are caused by an energetic event in Westerlund 2, presumably due 
to an anisotropic supernova explosion of one of the most massive member stars. While 
the origin of the TeV and GeV gamma-ray sources is uncertain, one may speculate 
that they are related to the same event via relativistic particle acceleration by strong 
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shock waves produced at the explosion or by remnant objects such as a pulsar wind 
nebula or microquasar. 

Key words: ISM: molecules, stars: individual (Westerlund 2), acceleration of 
particles, stars: supernovae: general 

1. Introduction 

Massive stars have a strong influence on the surrounding interstellar space via their stellar 
winds and ultraviolet radiation. Moreover, they end their lives in catastrophic supernova (SN) 
explosions, thus providing more energetic impact on the interstellar medium. Wd2 is a rich star 
cluster of age 2 — 3 Myrs (Piatti et al. 1998). CO observations have allowed a distance estimate 
to the source of 5.5 ± 1.5 kpc (Furukawa et al. 2009; see also Dame 2007). The total stellar 
mass in Wd2 is of the order of 45OOM0, including 12 O stars and 2 WR stars (Rauw et al. 2007). 
It is also associated with a HII region RCW49, a remarkable infrared nebula as revealed by 
Spitzer (Churchwell et al. 2004), and perhaps with the extended TeV gamma-ray source HESS 
J1023-575, discovered with the HESS telescope array in the direction of the cluster (Aharonian 
et al. 2007). Recently, low energy (MeV/GeV) gamma-rays have been reported from the same 
direction by the Fermi collaboration (Abdo et al. 2009). 

In this Letter we describe the discovery of a spectacular jet and arc of molecular gas 
detected with the NANTEN telescope in the ^^C0( J =1 — 0) mm- wave emission line survey 
(Mizuno & Fukui 2004), and discuss their possible links to the MeV/GeV and TeV gamma-ray 
sources. 

2. Results 

We show the two new CO features toward Wd2 in Figure 1, where the CO distribution 
in the velocity range ^ 20 — 30 km s~^ is presented in Galactic coordinates. One is a straight 

feature in the east at (/,6) = (284°. 2 - 284°.9, -1° .0 0°.4) (hereafter, the jet) and the other 

an arc-like feature in the west at (/,6) = (283°.9 - 284°.3, -0°.5 - 0°.0) (hereafter, the arc). We 
note that the arc has an intensity peak and a protrusion to the west of HESS J1023-575 at 
(/,6) = (284°.04-284°.09,-0°.31--0°.24), both of which lie along the jet axis, suggesting the 
existence of the diametrically opposite part of the jet. 

These two features show a remarkable correlation with Wd2 / HESS J1023-575. The 
jet is extended over 1° to the east of Wd2 / HESS J1023-575, tilted at an angle of - 39° to 
the galactic plane, and is very well aligned with Wd2 / HESS J1023-575 as well as a ridge of 
the radio continuum corresponding to the HII region, RCW 49(Whiteoak & Uchida 1997). A 
least-squares fitting to the jet in Figure 1 yields a regression, 

b + r .0 = (-2:50 ± 0:20) + (0.81 ± 0.01)(/ - 285°.0), (1) 

2 



where / and b are in degrees. 

This hne passes through the la error box of the peak position of HESS J1023-575, which 
is shifted from the center of the cluster by ~ 6 arcmin. The arc is apparently symmetric with 
respect to the jet axis and is well centered on HESS J1023-575. It appears somewhat rim- 
brightened along the outer boundary and shows a marked crescent shape, suggesting part of a 
swept-up shell. For convenience we define an axis which runs along the direction of the jet 
and passes through the peak position of HESS J1023-0575, and an axis T which is perpendicular 
to S. Here the origin of both T and S is = (284M9, -0°.39). 

The velocity distributions of the two features are shown in two overlaid velocity channel 
maps (Figure 2a) and in a position- velocity map along S (Figure 2b). 

The jet is generally quiescent, having a velocity width of a few km s~^ between S ^ 
0°.l — 0°.6 but show fairly large widths of 10 km s~^ at S ^0°.5 and 0°.7, if all the features are 
assumed to be physically associated. We tentatively divide the molecular gas towards the jet 
into the following four components; 

Jl: A narrow component at S ^0°.0 — 0°.6 located at ^ 29 km s~^ with a linewidth of 1 
— 2 km s~^, with little velocity gradient. This component has three local peaks at S' ~ 0°.2,0°.3, 
and 0°.5 with a winding shape towards ^ 0°.l — 0°.3. 

J2: A broad component at S ^ 0°.6 — 1°.0 located over a velocity range of 19 — 30 km 

J3: A fairly broad and localized feature at S' ^ 0°.5 and I^lsr 21 — 24 km s~^, towards 
one of the peaks of Jl. 

J4: The western component partially overlapping with the arc at S' ~ — 0°.2 and V^sr ^ 26 
km s"-*^. 

Although it is not conclusive, we suggest that the four components of the jet are physi- 
cally related as discussed later. The ^^CO linewidth of the arc is ^ 4 km s~^ centered at 26 km 
s~^, and we infer that the possible expansion velocity of the arc is not large. 

We shall here assume that the jet and arc are at ~ 5.4 kpc with a 30% error limit, and 
estimate the mass of each as follows, where an X factor of 2.0 x 10^^ cm~^ (K km s"-*^) is used 
to convert the ^^CO intensity to molecular mass (Bertsch et al. 1993): M(J1) ^ 7.9 x lO^M©, 
M(J2) - 2.2 X 10^ M0, M(J3) - 3.0 x 10^ M©, M(J4) - 1.0 x 10^ M©, and M(arc) - 2.0 x 10^ Mq, 
respectively. The accuracy in mass is limited to ±60% because of the uncertainty in distance. 

We show the HI distribution at I^lsr = 25.56 km s~^ with a width of 0.82 km s~^ in 
Figure 3 (McClure-GriflBths et al. 2005). The HI shows a clear sign of a hole towards Wd2 
and HESS J1023-575. Hereafter we shall call the HI feature surrounding the hole the HI shell. 
This shell is seen over a velocity range of 18 — 32 km s~^. We note that the velocity width 
of the hole is rather small, a few km s~^, suggesting that the possible expansion velocity is 
small at present. The HI hole has an elliptical shape with dimensions of ^ 29 pc by ~ 14 
pc, elongated along the direction of S. It is noteworthy that the HI shell exhibits an intensity 
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Fig. 1. The distribution of the ^^CO emission integrated over a velocity range from 24 — 28 km s~^ (the 
arc in the dashed region) and from 28 — 30 km s~^ (the jet). The orange contours correspond to the 
TeV gamma ray source, HESS J1023-575 (Aharonian et al. 2007). The red circle indicates the location of 
the gravity centre of the low energy (MeV/GeV) gamma-ray source reported by the Fermi collaboration 
(Abdo et al. 2009). The white contours are the radio continuum (Whiteoak & Uchida 1997). The cross 
indicates the position of Wd2. The axis S is defined along the direction of the jet with the position of 
HESS J1023-575. 

depression in its northern part, coincident with the molecular arc, suggesting that the arc and 
shell are physically related. This HI intensity depression is likely due to the conversion of HI 
into H2. The jet is coincident with an elongated spur of HI extended to the east over ^ 1°, 
from / ^ 284°. 3 to 285°. 0, while on the west there is little HI beyond the arc. The HI mass of 
the cloud shown in Figure 3 amounts to ^ 1.5 x IO^Mq in the velocity range 18 — 32 km s~^ 
for a conversion factor of 1.8 x 10-^^ cm~^ (K km s~-^)~-^ (Dickey et al. 1990). 

3. Discussion 

The CO arc and HI shell suggest a spherical shock either due to a SN or stellar winds 
of WR stars. The secure lower limit for the kinetic energy involved in the arc is estimated to 
be ^ 2.0 X 10^'' ergs for a molecular mass of ^ 2.0 x 10^ and an expansion velocity of 1 km 
s~^. A more likely estimate is perhaps ^ 1 x 10^^ ergs or more if we add the expansion energy 
of the HI gas to the expansion energy of the gas traced in CO. A SN with mechanical energy 
^ 10^^ ergs is more than sufficient to account for the energetics of the observed molecular and 
atomic gas, even when taking into account the conversion efficiency into kinetic energy of ^0.1 
(Thornton et al. 1998). There are many O stars in Wd2, whose collective mechanical wind 
energy may be large. For example 10 0-stars each with stellar-wind kinetic energy luminosities 
of 10^^ erg s~^ yield ~ 10^'' erg s~^. Over a timescale of ~ 10^ years the total kinetic power could 
reach ~ 3 x 10^^ ergs. With the same conversion efficiency, 0.1, a kinetic energy of ^ 3 x 10^^ 
ergs is available. 

The one-sided molecular jet can be explained in terms of in-situ conversion of HI into 
molecular gas by shock compression as a result of the jet propagating through the atomic 




S (degree) 

Fig. 2. Velocity distribution of the arc and the jet along the axis of the jet. (a) Integrated intensity map 
in ^^CO emission. The velocity ranges are 24 — 28 km s~^ for the arc and 28 — 30 km s~^ for the jet. 
Dashed contours show emission in the range 25 - 28 km s"^ and are every 1.6 K km s ^ from 2.4 K km 
s~^ (~ 3a). The blue mark indicates the position of the TeV gamma ray source, (b) Position- velocity 
diagram taken along S in figure 1. The horizontal axis indicates the coordinate of S with the position of 
the TeV gamma ray source as the origin. The definitions of the components Jl, J2 and J3 are described 
in the text. 

medium over a projected length of ^ 100 pc. This is consistent with the huge mass of the 
molecular jet — several times lO^M© — and the considerable amount of HI gas 5.2 x 10^ Mq 
) seen towards it. A lower limit for the time scale of the molecular jet is estimated to be ^ 10^~^ 
yrs from the timescale of shocked CO formation (Koyama & Inutsuka 2000; Bergin et al. 2004). 
There are presently two other known candidates for molecular jets driven by high-energy jets — 
SS433 and MJG348.5 (Yamamoto et al. 2008). The physical parameters of the two molecular 
jets are estimated as follows; full length ^ 500 pc, molecular mass ~ IO^Mq, and minimum 
kinetic energy ^ 10^^ ergs. For the Wd2 clouds, we roughly estimate such kinetic energy to be 
^ 2.0 X 10^^ ergs from the mass, ~ 2.2 x 10^ M©, and the velocity span, ^ 10 km s~^, of the 
dominant component J2, which is somewhat larger than the above two cases. We note that an 
increased velocity dispersion at the far tip of the jet, as seen in J2, is also observed in SS433 and 
MJG348.5. This is interpreted as a result of the stronger dynamical interaction towards the tip 
of the jet, where the deceleration of the high-energy jet becomes most significant (Yamamoto 
et al. 2008). 

We now consider the origin of the jet and arc. Presently, there are two scenarios in which 
highly collimated flows of sufficient energy and length may be formed: (i) a highly anisotropic 
supernova explosion; (ii) a high-energy accretion-powered jet from a compact object such as in 
a microquasar. Anisotropic energy output is thought to have occurred in several known SNe 
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including Cas A and W49B (Hwang et al. 2004; Miceli et al. 2008). The expected speed of 
the jets is > 1000 km s~^, suggesting a travel time of ^ 10^ yrs over ^ 100 pc if the jet axis 
is nearly perpendicular to the line of sight. Theoretical calculations indicate that a strongly 
magnetized and rotating neutron star formed in a SNe results in highly collimated jet (Burrows 
et al. 2007). Such explosions are able to release > 10^^ ergs in energy (Komissarov & Barkov 

2007) , and the compact remnant is expected to be observable as a magnetar. The frequency 
of known magnet ars is low, with only 12 magnet ars (Mereghetti 2008) amongst the 1500 or 
so known pulsars (Manchester 2004), while jet-like features in historical SNRs such as Cas A 
suggest that the anisotropic SNe may not be very rare. To summarize, the anisotropic SNe 
scenario may well explain the jet and arc as caused by a single collimated SNe. 

An alternative is that a conventional, isotropic SNe occurred in a binary system, leading 
to the formation of an accretion-powered jet. For example, the microquasar SS433, exhibits 
an X-ray jet of 150 pc and molecular jet of 400 pc in length (Kotani 1998; Yamamoto et al. 

2008) . In Wd2 the arc and shell might be formed by the SN explosion, and the molecular jet 
by a long-lived microquasar jet. The SS433 jet has a momentum flow rate of 1.5 x 10"'' Mq 
yr~^ and can supply 10^^ ergs in kinetic energy in 10^ yrs (Kotani 1998; Marshall et al. 2002), 
sufliciently large to supply the kinetic energy above. 

Concerning the origin of the gamma-ray emission, the intriguing alignment of the gravity 
centers of both the Fermi and HESS source positions with the molecular jet and arc suggests 
a common (initial) event, e.g., a supernova explosion. The link could be direct, via radiation 
of relic particles produced at the SN explosion, or via on-going acceleration of particles by a 
relativistic object such as a remnant of the explosion, e.g. supernova remnant shell, a pulsar 
wind nebula (PWN) or a microquasar jet. The flrst scenario favors a hadronic origin for the 
gamma-ray emission. Because of severe radiative (synchrotron and inverse Compton) losses, 
the relic ultrarelativistic electrons hardly could survive to the present epoch. Indeed, for any 
reasonable parameters characterizing the ambient medium, the lifetime of multi-TeV electrons 
is signiflcantly shorter than the age of the system. On the other hand, inverse Compton 
radiation of gamma-rays by continuously accelerated electrons, e.g., by a PWN, is a viable 
option. The gamma-ray luminosity is about 10^^ erg s~^. If the main target for IC gamma-rays 
is the 2.7 K CMBR, this mechanism would require acceleration of TeV electrons at a rate of 
g-i^ order of magnitude less, if the optical radiation of stars plays a dominant 

role. The latter case can be realized when acceleration and radiation of electrons take place 
not far from the stars. The extended character of the gamma-ray emission excludes a single 
star origin for the observed TeV emission. In this case the extension of the gamma-ray source 
is effectively determined by the volume occupied by stars, and consequently the size of the 
source is expected to be energy-independent. In contrast to this scenario, for a single PWN we 
expect an energy-dependent morphology caused by radiative energy losses which do not allow 
the highest energy electrons to propagate to large distances. Such an effect is observed in the 
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PWN HESS J1825-137 (Aharonian et al. 2006). In the case of a hadronic scenario we expect 
just the opposite dependence. Since generaUy the higher energy protons propagate faster than 
low energy protons, we might expect an increase of the angular size of the source with energy 
(Aharonian 2004). Another test to distinguish between the hadronic and electronic models 
is the ratio of GeV and TeV gamma-ray fluxes. In the case of the IC mechanism, we expect 
suppression of low energy gamma-rays and a rather low GeV/TeV ratio. It is expected that 
future higher quality GeV and TeV data should allow us to conduct detailed a quantitative 
study of the morphological and spectral characteristics of the radiation in different energy 
bands, which could lead to deflnite conclusions. Here we would prefer to limit the discussion 
by noting that the total energy required to explain the TeV gamma-ray data reported by the 
HESS collaboration for a source located at a distance of 5 kpc is about 10^^(n/lcm~^)~^ erg. 
The CO observations indicate a rather low density of gas in the region that coincides with the 
location of the TeV gamma-ray source (see above), most likely, n < 1 cm~^. Then the total 
budget in relativistic protons is required as large as 10^^ erg. This can be considered as an 
argument in favor of a very strong SN explosion with total energy exceeding 10^^ erg. 

In summary, our observations conducted with the NANTEN telescope in the ^^C0( J = 
1 — 0) mm- wave emission led to the discovery of a spectacular jet and arc of molecular gas 
toward the young star cluster Westerlund 2 which may be the result of a powerful supernova 
explosion. Another consequence of this explosion could be the GeV and TeV gamma-ray sources 
located in the same region as reported recently by the Fermi and HESS collaborations. Further 
detailed studies of spectral and morphological features of gamma-ray emission are requited to 
explore the links between these two phenomena. 

The original NANTEN telescope was operated based on a mutual agreement between 
Nagoya University and the Carnegie Institution of Washington. This work is flnancially sup- 
ported in part by a Grant-in- Aid for Scientiflc Research (KAKENHI) from the MEXT and 
from JSPS, in part, through the core-to-core program and by the donation from individuals 
and private companies. 
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